Decoherence in a quantum harmonic oscillator monitored by a Bose-Einstein 

condensate 
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We investigate the dynamics of a quantum oscillator, whose evolution is monitored by a Bose- 
Einstein condensate (BEC) trapped in a symmetric double well potential. It is demonstrated that the 
oscillator may experience various degrees of decoherence depending on the variable being measured 
and the state in which the BEC is prepared. These range from a 'coherent' regime in which only the 
variances of the oscillator position and momentum are affected by measurement, to a slow (power 
law) or rapid (Gaussian) decoherence of the mean values themselves. 
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In the past few years there has been much interest, 
both theoretical and experimental, in nano-mechanical 
oscillators whose quantum behaviour can be observed 
(measured) within the limits imposed by the uncertainty 
relations [J, Typically, various degrees of coherent 
control over such oscillators can be achieved by incorpo- 
rating them in hybrid devices involving superconducting 
microwave cavities Q , superconducting qubits [3, Q , sin- 
gle electron transistors [f| 0] and point contacts (PCs) 
[8j]. More recently, several schemes for coupling a quan- 
tum system to a Bose-Einstein condensate (BEC) have 
been proposed In the case of a measurement involv- 
ing a PC in a large bias regime, interaction between an 
oscillator and the electron current damps the latter leav- 
ing it in an equilibrium thermal state [lfjj . In this Let- 
ter we analyse a setup in which an oscillator is coupled 
to a BEC trapped in a symmetric double well potential 
rather than to a PC. With the atomic current dependent 
on the oscillator coordinate, the BEC is able to monitor 
the oscillator evolution, at the cost of introducing de- 
coherence to the oscillator dynamics. This decoherence 
is the main subject of this Letter. We will show that 
an oscillator monitored by a BEC does not, in general, 
undergo a quantum to classical transition [TjJ, LU| and 
may, in some cases, retain a degree of coherence. For re- 
cent relevant work on the types of decoherence possible 
in open systems we refer the reader to Ref. [l2|. 

We consider a system described by the Hamiltonian 
which is a _generalisation of the 'gatekeeper' model intro- 
duced in [13| (we put H = 1) 




FIG. 1: Double-well potential containing the N bosons that 
may tunnel from one well to the other. The flow of bosons is 
modulated by some function of the position of the oscillator. 
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c L (c R ) are the creation operators for the particles of the 
condensate in the left (right) reservoir and H osc is the 
oscillator Hamiltonian. The operator A osc represents the 
variable which controls the tunnelling rate between the 
two wells, so that its evolution can be monitored by ob- 
serving the atomic current or a change in the number of 
bosons in one of the reservoirs. Coupled in such a man- 
ner, the BEC shares with a conventional von Neumann 
meter [Tij the property that if prepared in a stationary 
state, it will remain in that state inducing an additional 
force on the measured oscillator [lfij]. With H con and 
Scon commuting, [H con , B con ] = 0, such states are easily 
found to be 
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|0) con , (n = 0,...,N), 
(3) 

where the vacuum |0) con corresponds to no bosons in the 
condensate. Assuming that the oscillator and the BEC 
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are prepared in a product state, p(0) = cr osc (0) ® p C on(0), 
and noting that P C on|</>n) = (-/V — 2n)|0„), we find the 
state of the monitored oscillator at a time t by tracing 
out the BEC degrees of freedom, 



Pn = Tr c 



[|<M(^n|Pco„(0)], 



where 
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(5) 



"H osc (e) = P osc + eA) SC , and e„ = <5fi(iV - 2n). 

Thus, CTosc(i) is an incoherent superposition of the 
states obtained by evolving ct osc (0) with the family of 
Hamiltonians "H osc (^n), n — 0,1, ...,N, weighted by the 
probabilities P n to find the BEC in the state \<j> n )- Ac- 
cordingly, at a time t, the expectation value of an oscilla- 
tor variable represented by an operator O osc is given by 
the sum 



(Oosc) 



P n Tr n 
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Following [13| we take the limit in which the number 
of atoms becomes large, while the coupling between the 
oscilllator and each individual atom is reduced, namely 



N 
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and choose = so as to exclude a constant back- 
ground current. The conditions Q ensure that a macro- 
scopic atomic current flows from the left to the right well, 
while the Rabi period of an individual atom tends to in- 
finity. Thus, the atoms are not going to return to their 
initial state in the foreseeable future, i.e., the BEC be- 
comes an irreversible meter with a large number of de- 
grees of freedom [T^]. Replacing sums by integrals, 
25QJ2 n ~^ I-oo we rewr ite © as 
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where Ef and \tp\) are the eigenvalues and eigenvectors 
respectively of the Hamiltonian 'Hosc(e) and P(e n ) = 
P n /(25£l). The long time behaviour of cr osc (t) now de- 
pends on the spectra E\. Indeed, for E\ — Ej ^ const(e) 
rapidly oscillating exponentials will cause the second 
term in Eq. {[8j) to vanish, so that a osc (t) [and with it the 
averages ©] will tend to stationary values as t — >• oo. 
Without such a cancellation, the oscillator will not be 
able to reach a steady state no matter how long one waits. 

Consider further a BEC initially localised in the left 
well, p con (0) = |0 O )(M where |0 O ) = (cl) N \0) con /VW., 
coupled to a harmonic oscillator of mass m and frequency 



Do, H osc — uj (a^a + |) = P 2 /2m + ^muj 2 X 2 , charged 
so as to affect the barrier between the two wells and the 
tunnelling rate. Both couplings linear (A osc ~ X) and 
quadratic (A osc ~ X 2 ) in the oscillator coordinate are 
possible [TjI]. For the probability weights P n in Eqs. © 
and ((6]) we have 



P n = N\/(2 N (N-n)lnl), 



(9) 



and applying the Stirling formula in the limit ([7]) yields 
P(e) = (2-KK 2 y 1 / 2 exp(-e 2 /2k 2 ). Rather than analyze 
the density matrix © it is convenient to consider the 
mean position and momentum of the monitored oscil- 
lator, together with their variances, thus choosing the 
operator O osc in Eq.© to be X, X 2 , P or P 2 . 

A. Coherent motion with 'breathing'. With A osc = 
a^+a = \ / 2mujoX = X/Xq the BEC monitors oscillator's 
position X, and we need to consider motion in a family 
of harmonic potentials shifted relative to the original one 

by Sx n e 



0,1, ...,N, 



V n (x) — mui (x + 5x n ) 2 /2 — mu; 2 .5x 2 l /2. 



(10) 



Thus, the energy differences in Eq.© are independent 
j), and no steady state can be 
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reached. We note further that since % OS c( e ) m Eq.© 
remains quadratic in both X and P, equations of mo- 
tion for the five operators X, P, X 2 , P 2 and XP + PX 
form a closed system which can be solved for each value 
e„. Averaging the results with the probabilities P„ (c.f. 
Eq.©) then yields 



d(X(t)) 
dt 

dim 

dt 



= -<P(t)) 
m 

= -mui 2 (X(t)) - y / 2moj o y^^ P n e n . 
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Probability distribution in Eq.© is symmetric about 
N/2 so that the sum in the last equation vanishes, and 
we find the mean values of both the coordinate and 
the momentum unchanged by the presence of the BEC, 
(X(t)) = (X(t)) free and (P(t)) = (P{t)) free , where the 
subscript 'free' refers to an oscillator uncoupled from the 
BEC. This does not, however, imply that the BEC has 
no effect on the dynamics of the oscillator. Indeed, cal- 
culating the variances we find 



(AX) 2 ree 



(AX) 2 = (X 2 ) - (X) 2 
(AP) 2 = (P 2 ) - (P) 2 = (AP) 2 f 
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where a x ee , U- M and a P ee V^o^. 



Thus, while (X(t)} and (P(t)} follow their unperturbed 
trajectories, the widths of the corresponding distribu- 
tions 'breath', first increasing and then decreasing again. 
Figure O shows the dynamics of the corresponding mean 
values and variances for an oscillator prepared in a coher- 
ent state (minimal Gaussian wavepacket) (a;|V , osc(0)) = 
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we find 




FIG. 2: Coherent motion (X is monitored) of a coherent 
initial oscillator state with (X(0)} /Xo = 0, <P(0)) /Po = 
2; X = (2mo) )" I/2 , Po = (mw /2) 1/2 . a) Mean po- 
sition (X(t)) /Xo (thick solid) vs. cjot. Also shown are 
[(X{t))±AX]/X (solid) and [(X{t))±AX free ]/X (dashed); 
b) Mean momentum (P(t)} /Po (thick solid) vs. u>ot- Also 
shown are [(P(t)} ± AP]/P (solid) and [(P(t)) ± AP free ]/P 
(dashed). Inset: (P(t)) /P vs. (X{t))/X . 



( mU j / % )V± e -™olx-(X(0)))V2 e i(P(0))x with 5 fi 2 iV/^ = 

25. Note that AX(t) recovers its original value after ev- 
ery period T = 2tt /ujq as the oscillator returns to its ini- 
tial state in each V n (x). The momentum variance AP(i) 
does so also after every half-period, when the shape of 
the original wavepacket is restored but the position of its 
centre is reflected with respect to the origin of each V n , 
i.e., when AX(t) reaches its maximum value. 



B. Gaussian decoherence. With A a 



X 2 /Xq the BEC monitors the square of the oscillator's 
position, X 2 , and we need to consider motion in a family 
of harmonic potentials with the same origin, but with 
different frequencies, cj„ = \/uJq + 4e„wo, n = 0, 1, N, 



V n (x) = m{uj 2 + Ae n uj )x 2 /2 = mw 2 X 2 /2. 



(12) 



Now E\ — Ej = (i — j) y^uj 2 + 4ew 7^ const(e), and we 
expect that in the irreversible limit ([7]) oscillator will un- 
dergo a relaxation to the steady state given by the first 
term in Eq.([5]). Solving equations of motion for each os- 
cillator frequency ui n and averaging with the probabilities 



(X(t))=^P„<j<X(0))cos(w n t) 

n 

{x 2 {t)) = Y. p i 



(P(0)) 
mco n 



sm(u: n t) 



((XP + PX)(0)) . . , 
v -sm(2w„t) 
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(P(t)) 



(1 - cos(2w„i)) + 
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(1 + cos (2uj n t)) 
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((XP + PX)(0)) 



-=Y J Pn cosKt) - (X(0)) sinKi) 



sin (2u> n t) 



+ ^Pl (1 + cos (2w n t)) + (1 - cos (2w n t)) 

where (Oosc(0)) = Tr osc [O osc cr osc (#)] is the expectation 
value of O osc in the initial oscillator state. We note that 
for e„ < — wo/4, oj n becomes imaginary as the interaction 
turns the oscillator potential into a parabolic repeller, 
leading to the break up of the system. Choosing the in- 
teraction to be small enough to neglect the possibility of 
breakup and taking the limit ([7]), we replace the sums in 
Eqs. (fT3"|) by integrals to obtain 



(X(t)) = (X(0))Re[/ (wo*)] 



TO) 



Im[/i(«ot)], (14) 



(P(t)> - (P(0))Re[/ (wot)]-mwo<A:(0))Iin[/_i(wot)] 
with 



//3=o,±i(t; cr) = 



1 
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/1+z 



(1 + z)/ 3 / 2 ' 

(15) 

where a = Ak/ujq and t = wot- For a weak coupling, 
a << 1, and times not exceeding l/wocx 2 , the exponent 
in (fT5")l can be expanded up to the second order in z. 
Replacing the lower limit of integration by —00 and eval- 
uating Gaussian integrals yields 

fpfc a) = (l + iT<j 2 /±)- l l 2 e- T ' 2 ° 2 l^ +2 ^^e^ + 0(a), 

(16) 

It is readily seen that, irrespective of the value of j3, 
\fp(ui{jt)\ decays on a time scale T = l/(uio<j) ~ ft -1 
so that for T « t ~ w /k 2 both the mean position 
(X(t)) and the mean momentum (P(i)) will have de- 
cayed to zero, regardless of their initial values, the de- 
cay being Gaussian in time. Figure [3] illustrates Gaus- 
sian decoherence for the same initial coherent state as 
in Fig. [U and for N = 100 and 5Q/u> = 0.0024. We 
note [and this is a general effect of a weak coupling, 
n « ujq, acting over a long time t » T , c.f. Eq. (fT3|) ] 
that in the final steady state of the oscillator the ini- 
tial (kinetic) energy is shared equally between its ki- 
netic and potential components, i.e., for t » T we have 

^ = ^ 2 { x 2 (t)) = iqgm. 
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FIG. 3: Gaussian decoherence (X 2 is monitored) of a coher- 
ent initial oscillator state with (X(0)) /X = 0, (P(0)) /P = 
2; X = (2mw )" 1/2 , Po = (mu /2) 1/2 . a) Mean po- 
sition (X(t)) /Xo (thick solid) vs. uiot. Also shown are 
[{X{t)) ± AX]/X (dashed); b) Mean momentum (P(t)) /P 
(thick solid) vs. w t. Also shown are [(P(t)) ± AP]/P 
(dashed). Inset: (P(t)} /P vs. (X(t)) /A given by Eqs. (fT3)l 
(thick solid) and p4!16p (dashed). 



C. Power law decoherence. Other types of decoher- 
ence are possible with different choices of the initial 
state of the BEC. For example, a power law decoher- 
ence can be achieved by replacing in Eq. (fl~5| the smooth 
Gaussian factor by a discontinuous one. Thus, choosing 



P n EE |(0„|0o)| 2 = e-^/E n >N/2- 

and zero otherwise and taking the limit 
with (X (0)) = yields 



for n > N/2 
3 for a state 



(X(t)) = (P(0)> 



asin(wot) + icos(wot) 
m(aw + l)(l + i 2 /a 2 )' 



(17) 



so that for t >> a, (X(t)) tends to zero as 1/t. 

In summary, in a hybrid setup involving an oscillator 
and a BEC in a symmetric double well potential we have 
an example of a quantum detector with a large number 
of degrees of freedom. In the irreversible limit the me- 
ter provides unidirectional macroscopic atomic current 
whose magnitude depends on the oscillator's position. 
Unlike in the case of a point contact, the measurement 
does not lead to universal damping of the oscillator and 
eventual thermalisation of its initial state. Rather, de- 
pending on the oscillator variable being monitored as well 
as on the initial state of the BEC, the oscillator may o 
may not undergo relaxation to a steady state and retain 
a degree of initial coherence. The absence of a quantum 
to classical transition predicted, for example, for an os- 
cillator coupled to a PC, is a consequence of the fact that 
a single energy level, rather than a broad energy band, 
is available for each tunnelling boson. 
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